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While 203 K high temperature superconductivity in H 3 S has been interpreted by BCS theory 
in the dirty limit here we focus on the effects of hydrogen zero-point-motion and the multiband 
electronic structure relevant for multigap superconductivity near Lifshitz transitions. We describe 
how the topology of the Fermi surfaces evolves with pressure giving different Lifshitz-transitions. 

A neck-disrupting Lifshitz-transition (type 2) occurs where the van Hove singularity, vHs, crosses 
the chemical potential at 210 GPa and new small 2D Fermi surface portions appear with slow 
Fermi velocity where the Migdal-approximation becomes questionable. We show that the neglected 
hydrogen zero-point motion ZPM, plays a key role at Lifshitz transitions. It induces an energy 
shift of about 600 meV of the vHs. The other Lifshitz-transition (of type 1) for the appearing of 
a new Fermi surface occurs at ISOGPa where new Fermi surfaces appear at the F point of the 
Brillouin zone here the Migdal-approximation breaks down and the zero-point-motion induces large 
fluctuations. The maximum Tc = 203A occurs at 160 GPa where EfIwq = 1 in the small Fermi 
surface pocket at F. A Feshbach-like resonance between a possible BEC-BGS condensate at P and 
the BGS condensate in different k-space spots is proposed. 

PACS numbers: 74.20.Pq,74.72.-h,74.25.Jb 


Introduction 

The recent discovery of superconductivity in pressur¬ 
ized sulfur hydride metal with critical temperature with 
Tc above 200K [Ei has provided experimental evidence 
that the coherent quantum macroscopic superconducting 
phase can occur at a temperature higher than the low¬ 
est temperature ever recorded on Earth (-89.2 °C) [1,0|- 
High temperature superconductivity has been found in 
other h;^rides @ like pressurized PH^ with Tc at about 
100 K 0, and it has been predicted in other hydrides 
like yttrium hydride [E- Superconductivity in pressur¬ 
ized hydrides was proposed by Ashchroft and his collab¬ 
orators Bi- The Universal Structure Predictor: Evo¬ 
lutionary Xtallography, USPEX, code [l^ available to¬ 
day has allowed to predict the structure and high tem- 
perature superconductivity in pressurized sulfur hydride 
prl - ll^ . The prediction of metallic H 3 S with the /m3m 
lattice symmetry [E [3 has now been confirmed by 
x-ray diffraction experiments Q above 120 GPa while 
at lower pressure different stoichiometry and structure 
are found [alii. The low mass of H atom is push¬ 
ing the sulfur-hydrogen T 2 U stretching mode and the Tiu 
phonons at P to high energy, so the energy cut-off for 
the pairing interaction ojq =150 ±50 meV . The super¬ 
conducting temperature was predicted by employing the 
Alien-Dynes modified McMillan formula [lll - H^ , and the 
Migdal Eliashberg formula [E-[E|- Most of these works 
have used the Eliashberg theory with BCS approxima¬ 
tions for isotropic pairing in a single band metal, assum¬ 
ing the dirty limit, which reduces a multiband to an ef¬ 


fective single band metal, and the Migdal approximation 
oJo/Ep « 1. More advanced theories have used density- 
functional theory including calculations of the effective 
Coulomb repulsion [23,[2lj. Within the Eliashberg theory 
the Migdal approximation assumes that the electronic 
and ionic degrees of freedom can be rigorously separated 
in agreement with the Born-Oppenheimer approxima¬ 
tion, which is valid for metals where the chemical poten¬ 
tial is far away from the band edges. The breakdown of 
Migdal approximation was observed in the multi-gap su¬ 
perconductor AlxMgi-j;B 2 pE - [E| where the band edge 
of the a band fluctuates across the chemical potential due 
to zero point motion [ 2 ^ The breakdown of the Migdal 
approximation in a multigap superconductor is relevant 
since it requires the correction to the chemical potential 
induced by pairing below the critical temperature [E| ig¬ 
nored in the standard Eliashberg theory. Moreover where 
the Migdal approximation breaks down in a multigap 
superconductor it enetrs in a unconventional supercon¬ 
ducting phase where the coexisting multiple condensates 
could be either in the Bose Einstein condensate (BEC) 
regime or in the BEC-BCS crossover regime. In these 
complex multi gaps superconductors where BEC, BEC- 
BCS and BCS condensates coexist the excha.nge inter¬ 
action between different condensates [E, [E - [E | which 
is neglected in the BCS approximations, but could be¬ 
come a relevant term both to increase or to stabilize high 
temperature superconductivity. This quantum term is a 
contact interaction, given by the quantum overlap be¬ 
tween the condensates, which increases the condensation 
energy [E HEHIl and the critical temperature via the 
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shape resonance, analog to the Fano Feshbach resonance 
in ultracold gases. It has been found that the critical tem¬ 
perature shows the maximum amplification where one of 
the condensates is in the BCS-BEC crossover, which can 
occurs on the verge of a Lifshitz transition [3, 0, 

At the Lifshitz transition, a change of the topology of the 
Fermi surface is induced by pressure or doping and it has 
been shown to control high temperature superconductiv¬ 
ity in iron pnictides [30|, |3l| . 

Since H^S is a multiband metal and Lifshitz transitions 
could occur by increasing pressure ii it has been pro¬ 
posed that it is a multigap superconductior near Lifshitz 
transitions (?7l - [^ . where also multi scale phase separa¬ 
tion at a the Lifshitz transition could a ppear, s imilar 
to what has been observed in the cuprates [34 - 1^ . 


The coupling between the electronic and the atomic 
lattice degrees of freedom in sulfur hydrides at zero tem¬ 
perature has been neglected in previous calculations of 
the electronic structure assuming a very large Fermi en¬ 
ergy. On the contrary in the case of band edges close to 
the chemical potential and near Lifshitz electronic topo¬ 
logical transitions the zero point motion (ZPM) cannot 
be neglected. It is known that the ZPM modifies the band 
structure itself [ssl - l^ . with corrections of the band gap 
energy which can be larger than those induced by corre¬ 
lation. Moreover, the lattice disorder from ZPM will also 
cause band broadenings, as it has been demonstrated in 
early investigations . Such effects have been shown to 
be important for many properties in several different ma¬ 
terials, even if their atomic masses are larger [42l - l43| and 
the lattice disorder can perturb spin waves and phonons 
( 4 ^ . The zero point motion is expected to be large in H 3 S 
because of the small mass of H atoms, high frequency S- 
H stretching modes, and the double well potential for 
hydrogen in H 3 S Q which is similar to the well studied 
cases of ice and biological macromolecules. 

In this work we discuss the effect of lattice compres¬ 
sion and zero point motion on the electronic structure 
of H 3 S in the pressure range above 120 GPa where 
the metallic Imim phase is stabilized by the ZPM of 
hydrogen atoms. We first focus on the effect of the 
lattice compression on the large Fermi surface, identified 
with No.2 in ref. Q and on the small hole Fermi 
surface pockets at F identified with No. 4 and No.5 in 
reference Q. In particular we discuss the character of 
van Hove singularity (vHs) in the large Fermi surface, 
giving a peak in the density of states, and how its moves 
toward the chemical potential by increasing pressure. 
Above the critical pressure 210 GPa the vHs crosses 
the chemical potential and we find that tubular Fermi 
surface portions with two-dimensional character in the 
k-space appear, indicating a 3D-2D topological Lifshitz 
transition of type 2. Second, we study the effect of the 
zero point motion on the van Hove singularity and we 
find a very large energy shift of the order of 600 meV 
of the vHs due to ZPM. Third in the pressure range 
120-180 GPa the amplitude of ZPM of hydrogen-sulfur 
stretching mode is shown to be larger than the S-H bond 


splitting in the R3m structure therefore it stabilizes 
the symmetric /m3m structure. In this pressure range 
we have studied the Lifshitz transitions of type 1, for 
the appearing of new small Fermi surface pocket at 
F at 130 GPa. Finally we compare our results with 
recent experiments showing that the isotope coefficient 
for the superconducting critical temperature diverges 
and the critical temperature goes toward zero at 130 
GPa in agreement with predictions of the BPV theory 
[27l - l3ll| for a Lifshitz transition of type 1. Moreover the 
maximum critical temperature, 203 K, appears where 
the Fermi energy in the small Fermi hole pocket at F is 
of the order of the pairing interaction wq as predicted in 
ref. [IMIl- 


The Band Structure 

The high pressure phase of metallic H 3 S has the cubic 
Space Group: 229 with ImZm lattice symmetry. The 
Im3m lattice structure can be described by the small 
Body Gentered Gubic (bcc) unit cell, which has been used 
by all previous calculations providing the electronic band 
dispersion in the bcc Brillouin zone (bcc BZ). 

Here, we use also an alternative simple cubic unit cell 
made by 8 atoms per unit cell with a simple cubic Bril¬ 
louin zone (sc BZ), which permits an easier comparison 
with a traditional group of superconductors, namely the 
A15 compounds. In fact A15 compounds have a lattice 
structure with the cubic Space Group: 223 belonging to 
the same ditesseral central class, or galena type, of cu¬ 
bic space groups with the same Hermann-Mauguin point 
group mZm. 

We have obtained the band dispersion in the large sim¬ 
ple cubic (sc) Brillouin zone which grabs more details of 
the electronic band dispersion in the complex bitruncated 
cubic honeycomb lattice. Here S sites form a bcc lattice 
(exactly as Si in the A15 compound V 3 Si) with linear 
chains of H on the limits of the sc cell, similar to how the 
transition metals form linear chains in A15 superconduc¬ 
tors. These results allow a more easy comparison of the 
electronic bands between the two type of materials. 

We show the electronic self-consistent paramagnetic 
bands for the bcc BZ in Fig. 1 and in Fig. 2 for the sc 
BZ. In this last picture the simple cubic unit cell contains 
8 sites totally. The calculations have been performed by 
the linear muffin-tin orbital (LMTO) method and the 
local spin-density approximation (LSDA) in the frame of 
standard methods [TTH^ . The basis set goes up through 
^=2 for S and £ = 1 for H. We have found that the S- 
2 p core levels are always far below the valence band re¬ 
gion at these pressures (about 12.2 Ry below Ep). The 
Wigner-Seitz (WS) radii are 0.38a for S and 0.278a for 
H. The k-point mesh corresponds to 11 points between F 
and X, or 1331 points totally in 1/8 of the BZ, or with 
finer k-point mesh for plots of the bands along symme¬ 
try lines. One spin-polarized calculation is made for the 
largest volume starting from an imposed ferro-magnetic 
(FM) configuration. All local FM moments converge to 
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FIG. 1: (Color online) The band structure for H 3 S at sym¬ 
metry points for the small bcc unit cell with a=5.6 a.u., cor¬ 
responding to P = 210GPa in the bcc Brillouin zone (bcc 
BZ) 
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FIG. 3: (Color online) The total DOS for H 3 S with variable 
lattice constants between 5.6 and 6.2 a.u. 

the largest lattice constants, see Fig. [3l 

An approximate value for the pressure, P, calculated 
as a surface integral using the virial theorem [5Ij , is useful 
to get insight to the relative contribution from different 
atoms. These total pressures (P) amount to 0.98 and -0.3 
Mbar at the two extreme lattice constants. The partial 
S- and H-pressures increases much more on H than on S 
when the lattice constant is reduced. This together with 
the charge transfers, indicate that the S sublattice is more 
compressible than that of H. A charge transfer from S to 
H at high pressure will enforce the hardening of H. More 
precise values of P, shown in Table HI are obtained from 
the volume derivative of the calculated total energies, and 
we obtain 180 GPa for a=5.5 a.u. and 73 GPa for a=6.1 
a.u., which agree well with ref. [l^ . 

TABLE I: Lattice constants a, DOS at E_f (in units of 
states/eV/Gell), pressure P, the Hopfield parameter NP in 
eV/"jP for S and H, A, and Tc estimated from the McMillan 
formula. The phonon moments are 600 and 1800 K for S and 
H, respectively, and uiog is 1500 K. 


FIG. 2: (Color online) The band structure for H 3 S at symme¬ 
try points for the simple cubic large unit cell, in the simple 
cubic BZ at a pressure of 210GPa 


zero, which shows that FM and FM spin fluctuations are 
unlikely. The low-lying s-band on S is very similar to the 
Si-s band in the A15-compound VaSi [^, but the high 
pressure in H 3 S makes it to overlap with the S-p band. 
In contrast, the Si-p band in VaSi is separated from the 
Si-s band. The separation is never complete in Has', but 
it is visible as a dip in the DOS at 13 eV below Ep for 


a (a.u. 

) N{Ef) 

P (GPa) 

AP(S) NIPH) 

A 

Tc (K) 

5.4 

1.00 

300 

5.6 

4.3 

1.0 

145 

5.6 

0.99 

215 

5.5 

4.3 

1.0 

143 

5.8 

0.90 

150 

4.5 

3.4 

0.81 

109 

6.0 

0.88 

95 

4.1 

3.1 

0.73 

93 

6.2 

0.90 

50 

3.9 

2.8 

0.70 

73 


The electronic numerator NP of the electron-phonon 
coupling constant A = NP/K (where the force constant 
K = Mp is taken from experiment) is calculated in 
the Rigid Muffin-Tin Approximation (RMTA) [ 5 ^ [s^. 
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The matrix elements for d — /-scattering in S and p — d 
scattering in H are missing, which leads to some under¬ 
estimation of the total A. Note also that LMTO uses 
overlapping Wigner-Seitz (WS) spheres with no contri¬ 
bution to P from the interstitial region. The values of 
NP on H and S are reasonably large in comparison to 
the values in transition metals. The low local DOS is 
compensated by larger P and the scattering to p-states, 
while in the A15’s large scattering arises mainly from d- 
states. Our N{Ef), shown in Table HI are in agreement 
with the work of Papaconstantopoulos et al. (l6l | , which 
quite unusually are largest at low volume and large P. 
This is because of the band edges that cross Ep and be¬ 
come occupied at high P. But for even larger volume 
N{Ef) behaves normally again and it increases when P 
goes down, as shown by Papaconstantopoulos el al. [l^ . 
The matrix elements P show strong P-dependences, and 
NP (the so-called Hopfield parameter) increases steadily 
with P in agreement with ref. M- The absolute value 
of A is smaller than in ref. [l^. This can partly be due 
to the smaller basis in our case, and also because of the 
use of WS-spheres in LMTO instead of non-overlapping 
MT-sphere geometry in LAPW [l^ . 


In the estimation of P from the McMillan equation 
another uncertainty concerns the Coulomb repulsion p*. 
Many Tc-calculations use /r*=0.11-0.13, but theories for 
calculation of p* are approximate or unreliable . Re¬ 
tardation can make p* large and screening makes it small 
or even negative depending on the band width and the 
phonon frequency ^ . Here we use an empirical formula 
for pL* proposed by Bennemann-Garland |55l |. which leads 
to a small value of the order 0.03. Thus, we caution that 
even if our calculated P from the McMillan formula will 
be large it is very approximate, similar to what has been 
concluded in the other works. Our P's are of the order 
145 and 75 K between the two extreme lattice constants 
when the total coupling strength goes from 1.0 and 0.7. 
These P’s are reduced to 95 and 40 K between the ex¬ 
treme values of a if p* is 0.13. The different contributions 
to the total A from each S and H atoms are comparable. 
Moreover, the total A itself is not unusually large. Us¬ 
ing the McMillan formula, valid only for a single effective 
band BCS system in the dirty limit, the high P is mainly 
because of the large phonon frequency pre-factor {uiog) 
in the equation. 


In order to investigate the pressure effects on the elec¬ 
tronic structure we show in Fig. |3]the different total DOS 
for different lattice constant a of the perovskite structure 
changing from 6.2 to 5.6 a.u. The total DOS at Ep per 
S-atom is about 3.5 {Ry)~^ and 1 {Ry)~^ per H-site, 
compared to the order 20 per V in V 3 Si or in elementary 
V and Nb. This is not surprising because of the large 
width of the wide band in H 3 S having its bottom at 2 
Ry below the chemical potential, see Fig. [31 while for 
increasing number of (d-) electrons in transition metal 
A15 compounds the total band width is more like 2/3 of 
aRy [13. 

The charge within the H WS-sphere increases from 


1.3 to 1.4 el./H when the lattice constant decreases from 
6.2 to 5.6 a.u., and the H-s charge goes from 0.95 to 
1.0. This fact justifies somewhat the use of LSDA for 
H even though atomic H with exactly 1 electron is best 
described by the Hartree potential only. The results 
show a very large effect of the pressure on the lowest 
dispersive bands with H s character, and the energy shift 
of the narrow peak of DOS at the chemical potential 
due to the van Hove singularity. The shift is small in 
the figure because of the large energy scale. 

Hydrogen Zero Point Motion Effects 

Usually the atomic velocities (vi) from the vibrations 
are much slower than the electronic velocities, and the 
electronic structure can relax adiabatically at all times. 
Therefore, the electronic structure calculations usually 
can neglect Ui, as in ’frozen-phonon’ calculations. How¬ 
ever it is well-known that the lattice of real materials 
becomes distorted at large T (’thermal disorder’) and 
some distortion remains at r =0 due to the ’zero point 
motion’, ZPM. Thermal disorder and ZPM will modify 
and broaden the bands compared to the case of a perfect 
lattice because the potentials at different sites are not ex¬ 
actly the same. Each atomic (i) lattice position deviates 
from its average position by u(t,T), because of the exci¬ 
tation of phonons The time average oi u(T) for harmonic 
uncorrelated vibrations are well represented by a Gaus¬ 
sian distribution function with width (FWHM) < u >, 
which tends to ^Hlo/K at low T (ZPM) and to ikpT/K 
at high T, see ref. [s^. 

The Debye temperature for H phonons (~ ISOOAT) in 
sulfur hydride is much higher than P and the disorder 
amplitude from ZPM is almost the same as at T = P- 
The maximum superconducting gap (2A) would be ^ 65 
meV at T—0 in the large Fermi surface, i.e., the band 
No.2 using the notations to identifythe 5 bands crossing 
the Fermi level introduced in ref. Q- With the parame¬ 
ters as in the calculation of A, and A = U, the amplitude 
for the potential modulation for phonons [s^l, one can 
estimate that ua, the displacement of phonons that lead 
to superconductivity should be ^0.02 A, i.e., less than 
the amplitude of ZPM for H, or comparable to the ZPM 
of S. 

We have used a first-principle DFT approach based on 
a large supercell to calculate the effect of ZPM on the 
van Hove singularity. This method it is more suitable for 
complicated 3-D materials such as H 3 S (having a unit 
cell of 64 atoms) compared with alternative perturbative 
approaches for evaluation of the band energy changes as 
function of displacement (u) based on the Allen-Heine- 
Cardona (AHC) method, which is suitable for simple sys¬ 
tems like carbon nanotubes [H, [s^- While AHG type 
approach may seem more sophisticated, when it is ap¬ 
plied to complex 3-D materials it is more important to 
consider good statistics and rely on experimental infor¬ 
mation on force constants The important point 

of our approach is that a large supercell is used in order 
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FIG. 4: (Color online) The DOS for HgS for cubic 64-site 
supercells. The (blue) thin line is the DOS for the perfectly 
ordered supercell. The (red) heavy line shows the DOS for 
the disordered supercell with us=0.01a and Mif=0.033a Zero 
Point Motion 


to have a good statistics for the individual u of each atom 
so that the bands energy changes as function of increas¬ 
ing disorder can be determined properly. The average 
atomic displacements and their 3-D distribution are de¬ 
termined from the phonon spectrum and atomic masses 
and it does not make sense to attempt an ab-initio cal¬ 
culation of u using frozen zone-boundary phonons, and 
tight-binding method. 

Calculations for disorder within supercells with 64 
atoms for FeSi [i^, or even less (48 atoms) for purple 
bronze [i^, have shown that different generations of in¬ 
ternal disorder in the cells produce similar results (as long 
as the disorder amplitude < u > are the same). Sym¬ 
metry makes the bands identical in different irreducible 
Brillouin zones (IBZ) for ordered supercells, and these 
bands have their exact correspondence in the bands of 
the small cell. But this symmetry is lost if the atoms are 
disordered, and the bands have to be determined in half 
of the BZ. 

With a force constant K = Moj'^ oil eV/A^ we obtain 
an average amplitude < u > of the order 0.15 Afor ZPM. 
This is close to 10 percent of the H-H distance, which 
according to the Lindemann criterion suggests that the 
H sublattice is near melting [H^. The < u >-amplitude 
for the S sublattice is normal, because of its large mass, 
and it is probable that the rigid S-lattice is important for 
the stability of the structure in which the H-atoms are 
rather loosely attached to their ideal positions. 

The ZPM generates different shifts of the Madelung 
potential at different sites. This modifies the band en¬ 
ergies (Ae) and leads to ’fluctuations’ (there is a spread 


of the eigenvalues in different IBZ) of Ae, in particular 
at the band edges. The low-T energy band fluctuations 
in materials with narrower band widths are known to be 
about 20 meV for u in the range 0.03-0.04 A[^-(45l|. 
From an extrapolation of these values to the conditions 
in H 3 S (larger < u >) we estimate that the band energy 
fluctuation can be larger than 150 meV for H-bands. 

The band No.2 of H 3 S in Fig. [3] is very wide, about 
2 Ry for the high lying valence bands, and the S-s band 
overlaps with the S-p band. This makes the band dis¬ 
persion and Fermi velocities high in large portions of the 
k space, where the effects of energy band broadening in 
these points of the k-space is less important. However in 
portions of the K-space around the F-X-M path in the sc 
BZ, where the van-Hove singularity crosses the chemical 
potential the Fermi velocity is getting small and strong 
dynamical fluctuations controlled by the zero point lat¬ 
tice fluctuations are expected to be relevant. 

The electronic calculations for the ZPM in H 3 S at 
P=210 GPa, a =5.6 a.u., have been carried out using 
a large supercell in which the lattice is disordered. Each 
atom is assigned randomized displacements, Ux,Uy,Uz in 
such a way that the distribution of all displacement am¬ 
plitudes (I u I) has a bell shaped (Gaussian) distribution 
with the FWHM width equal to the averaged displace¬ 
ment amplitude (< u >) at was described before and in 
ref. [ 4 ^ 1 . In Fig. |T] we show the energy renormaliza¬ 
tion of the DOS due to the calculated zero point motion 
(ZPM) . Here we consider a 2x2x2 extension of the cubic 
unit cell with 64 atoms totally which permits to calculate 
u from 192 displacement vectors, which is a reasonably 
good statistics for calculating < u >. As was discussed 
above, it is sufficient to do the calculation for one disor¬ 
dered configuration when the supercells are large with at 
least 48 atoms. Because of the large mass difference be¬ 
tween S and H we here allow larger < u > for H (uh) than 
for S {us)- For small < u > it can be shown that cor¬ 
relation of vibrational movements and anharmonic terms 
are small But uh is large, and in the generation 

of a disordered configuration for M^f/o=0.05, as for the 
expected ZPM, several pairs of H come too close to each 
other. Therefore, in order to avoid large anharmonic ef¬ 
fects at this stage, we calculate the electronic structure 
for a supercell with Ms/a=0.01 (which is close to the ex¬ 
pected ZPM for S) and M///a=0.033 (which is 2/3 of the 
expected ZPM for H). 

The resulting DOS for the disordered lattice is shown 
in Fig. [4] The results show a strong effect of ZPM on 
the energy of the narrow peak of the DOS just below the 
chemical potential. This peak has a large fraction of H s- 
and p-states, while the peak further down (3-4 eV below 
Ep) has less H-character and is less affected by disor¬ 
der. The width of the DOS peak with large character of 
H-states gets wider in the disordered case, which is ex¬ 
pected, but we see also a large energy renormalization of 
this narrow peak. This peak which was mostly below the 
chemical potential in the ordered lattice is now pushed 
up above the chemical potential because of the lattice 
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FIG. 6 : (Color online) The partial DOS with different atomic 
and orbital symmetry for H 3 S at the lattice constant 5.4 a.u.., 
the (red) thin lines are for S and the (blue) heavy lines for H. 


FIG. 5: (Color online) The total DOS functions in Fig. [3] on 
a finer energy scale near the chemical potential as a function 
of the lattice constant. Each DOS curve is separated by 0.5 
{cell eV')”^ units for the sake of visibility. 

ZPM. The shift at the top of the peak is of the order of 
600 meV (see Fig. [d]) and the peak becomes broader in 
the disordered case. 

Our simple approach to treat the ZPM has been 
based on the assumption that the frozen disorder in 
our supercell calculation is representative for ZPM. If 
so we may ask how much ZPM would change in the 
RMTA value of A and the McMillan estimate of Tc- 
The total DOS at Ep is almost the same (7.3 and 
7.4 states/eV/cell for ordered and disordered case, 
respectively) as can be seen in Fig. |4l The A's and T'^s 
are calculated to be 0.88 and 0.86, and 122K and 118Ar, 
for ordered and disordered cases, respectively. Thus 
there is a small reduction of A and from disorder 
even if the two N{Ep) are comparable. However, one 
can note that the peak in the DOS has moved from 
being below the chemical potential in the ordered case 
to be above when ZPM is taken into account indicating 
a 600 meV energy shift of the van Hove singularity with 
zero point lattice fluctuations involving H atoms. In the 
remaining sections we discuss some details about the 
Lifshitz transitions as a function of pressure. 

Lifshitz Transitions as a function of pressure 

The energy shift of the van Hove singularity (vHz) can 
be followed by looking at the shift of the narrow peak 
in the total DOS near the chemical potential at different 
lattice parameters a which is shown in Fig. [S] 

This peak is mostly due sulfur orbital contributions 
as shown by the partial DOS functions of the vHs for 
the case of lattice parameter a = 5.6a.u. shown in Fig. 



FIG. 7: (Golor online) The partial DOS for H 3 S at the fixed 
lattice constant 5.4 a.u. from different bands crossing the 
chemical potential in the simple cubic Brillouin Zone. 

El The S-p orbital and S-d orbital contribute to about 
30 and 16 percent of the total DOS respectively, there¬ 
fore the vHs is mainly due sulfur atoms. Both the to¬ 
tal charge and the ^-character at on H are mostly 
s (60-65 percent), and the hybridization between the s- 
electron and states on the S atoms away from the chains 
is large. This large on-site hybridization is favorable for 
large dipole matrix-element contributions to the electron- 
phonon coupling. Fig. [7] shows the band decomposition 
of the total DOS in the cubic 8-site unit cell in the sc 
BZ into 4 bands classified as No.8, No 9, No. 10, No. 
11. The vHs and the peak in the DOS is due only to 
band 10 which gives the largest Fermi surface and it is 
due to the flat dispersing bands with low Fermi veloc¬ 
ity near X and M points giving a van Hove singularity. 
Fig. [5]shows the pressure dependence of the vHs in band 
10 which goes through the chemical potential at about 
210 GPa where a is about 5.6 a.u.. This result shows 
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FIG. 9: (Color online) The band structure for H 3 S along X — 
M — F — R for a for the simple cubic double cell (sc BZ) with 
the lattice parameter changing between 5.4 and 6.2 a.u. 


FIG. 8 : The evolution of the partial DOS from band 8-11 
in the simple cubic Brillouin Zone as function of the lattice 
constant. Only bands 9 and 10, shown by thin and heavy 
continues lines respectively, contribute much to the total DOS 
at the chemical potential Ep. 


that the van Hove singularity approaches the chemical 
potential at 210 GPa but it remains near the chemical 
potential in the energy range of the energy cut off of the 
pairing interaction in the pressure range showing high 
temperature superconductivity. Where the vHs crosses 
the chemical potential a Neck-Disrupting Lifshitz Tran¬ 
sition, of type 2, occurs. Here the topology of the large 
Fermi surface changes because of the appearing of small 
pieces of tubular 2D Fermi surfaces connecting the large 
petals as discussed in the reference 0. In these tubular 
portions the Fermi velocity is small therefore the Migdal 
approximation breaks down. On the contrary in the large 
petals the Fermi energy is much larger and the Migdal 
approximation is valid. 

In order to associate the crossing of the chemical po¬ 
tential of the narrow peak in the DOS with the Lifshitz 
transitions on the topology of the Fermi surfaces we have 
plotted the electronic bands in a narrow energy range 
near the chemical potential for both the sc BZ and bcc 
BZ in Fig. [9] and Fig. [101 These band plots show that 
several Lifshitz transitions appear for increasing pressure. 
From Fig. |9]it is seen that a local band maximum crosses 
Ep in the sc BZ at about 2/3 of the F — M distance when 
a is 5.8. a.u. (see Figs. [Hand Figs. [5]) and the same 
band crosses the chemical potential between N and H 
in Fig. [ 10 ] in the bands for the bcc BZ. This Lifshitz 
transition occurs at P = 210GPa. The energy differ¬ 
ence E 2 between the chemical potential and this local 



FIG. 10: (Color online) The variation of the band structure 
near the chemical potential plotted for the small bcc unit cell 
of one formula unit of H 3 S (bcc BZ) as a function of the lattice 
constants between 5.4 and 6.2 a.u.. The colors of the bands 
correspond with different lattice parameters as in Fig. 9. The 
labels of symmetry points are as for bcc BZ. 


band maximum which is associated with the vHs goes 
from -200 to -1-100 meV when a decreases from 6.2 to 5.6 
a.u.. This gives a neck disrupting Lifshitz transition in 
the Fermi surface at 210GPa where the neck disappears 
at low pressure in the N-H direction [ 2 ^ . The Fermi sur¬ 
face neck appears at the point where the narrow peaks in 
the DOS crosses the chemical potential. Another band is 
seen to be approaching the chemical potential in Fig. 
between X and M when the pressure goes up. However, 
this potential band crossing (which is not seen on a sym- 
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FIG. 11: (Color online) The upper panel a shows the S- 
H bond length in the pressure range between 120 and 180 
GPa where without the hydrogen ZPM, the R3m structure 
was expected to be stable with the S-H bond splitting into a 
long (blue squares) and short (violet squares) sulfur-hydrogen 
bond. The amplitude of the calculated ZPM of the S-H bond 
is indicated by the red error bars Bi- The S-H amplitude 
of the zero point motion, ZPM, is larger than the S-H split¬ 
ting in the range 130-180 GPa therefore in this pressure range 
the ZPM stabilizes the /m3m structure B in agreement with 
experiments B|. The panel h shows the Fermi energy Ef4 in 
the small hole Fermi pocket at F as a function pressure. The 
top of this band crosses the chemical potential at 130 GPa 
giving the Lifshitz transition of type 1 for the appearing of a 
new Fermi surface. The position of the vHs E 2 remains below 
the chemical potential but it remains in the energy range of 
pairing interaction. Panel c shows the variation of the ex¬ 
perimental isotope coefficient calculated from data in ref. B 
which shows a divergence from 0.3 at 180 GPa to 1.5 at 135 
GPa, which is not predicted by the BCS theory. The critical 
temperature decreases toward zero with a decrease of about 
60 K in a range of 30 GPa, beteen 160 GPa and 130 GPa 
which is not predicted by the BGS theory. Both phenomena 
are predicted by the general theory of multigap superconduc¬ 
tivity near a Lifshitz transition [3, Q . 


metry line in the band plots for the small bcc cell) will 
not reach Ep to make a FS pocket unless P is increased 
even more. 

In the pressure range 120 < P < 160 GPa, where the 
onset of superconductivity occurs followed by the rapid 
increase of the critical temperature there are Lifshitz 
transitions, of type 1, for the appearing of new Fermi 
surface spots at F. There are 3 bands pushed up by in¬ 
creasing pressure which cross the chemical potential and 
three small closed Fermi surfaces appear. In fact Fig. |9l 
Fig. [TUI show that the tops of these bands at F are all 
above the chemical potential at the highest P while the 
tops of these 3 bands bands are below Ep at low P. 

In Fig. 11 we have plotted in panel a the S-H bond 
length as a function of pressure and the amplitude of the 
calculated spread of this bond length due to the hydro¬ 
gen zero point motion indicated by the red error bars. 
The USPEX theory [I^, [i3| predicts the second order 
phase transition from ImSm to R3m structure at 180 
GPa where the sulfur atoms remain in the same sites of 
the bcc unit cell while hydrogen ions are frozen in one of 
the two minima of their double well potential due to the 
hydrogen bond like in ice structural transitions. On the 
contrary this transition is forbidden in the pressure range 
120-180 GPa since the quantum zero point motion ampli¬ 
tude which is larger than the difference between the short 
and long hydrogen bonds expected in the R3m structure, 
Therefore the ZPM stabilizes the Im3m structure in the 
pressure ran ge 130-180 GPa Q in agreement with recent 
experiments [ 2 |. 

The panel b in Fig. 11 shows the Fermi energy in the 
small hole Fermi pocket at F as a function pressure. The 
top of this band crosses the chemical potential for a pres¬ 
sure larger than 130 GPa. At this pressure the Lifshitz 
transition of type 1 for the appearing of a new Fermi sur¬ 
face occurs. In fact above 130 GPa a new small Fermi 
pocket appears at the F point of the Brillouin zone. The 
Fermi energy Ep/i remains smaller than the energy cut off 
of the pairing interaction below 160 GPa. Therefore In 
the pressure range between 130 and 160 GPa the Migdal 
approximation breaks down for electron pairing in the 
small Fermi pocket at the F point of the Brillouin Zone. 

Finally panel c of Fig, 11 shows the variation of the 
experimental isotope coefficient in this pressure range 
calculated by recent data B which shows a regular 
increase with decreasing pressure from 0.3 at 180 GPa to 
1.5 at 135 GPa. The divergence of the isotope coefficient 
approaching the Lifshitz transition at 130 GPa is not 
predicted by the BGS theory using standard Midgal 
approximation and an effective single band but it is 
predicted in the frame of general theory of multigap 
superconductiivity near a Lifshitz transition Q This is 
supported by the decrease of the critical temperature of 
about 40 K in and of about 60 K in £>38 while in 
the BGS calculations the variation of the critical tem¬ 
perature over this 30 GPa range is predicted to be of the 
order of 10 K. In fact in the multigaps superconductors 
at a Lifshitz transition the critical temperature goes 
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toward zero (like at the Fano-Feshbach anti-resonance) 
at the appearing of a new nth Fermi surface where 
electrons in the nth band have zero energy Epn =0 at 
the band edge forming a BEC condensate, while the 
critical temperature has a maximum (Fano-Feshbach 
resonance) where the electrons have a Fermi energy of 
the order of the pairing energy forming a condensate in 
the BEC-BCS crossover Q. 

Conclusions 

In H 3 S the onset and the maximum superconducting 
critical temperature, 203 K, are controlled by pressure, 
like in cuprates where the onset and the maximum value 
of the critical temperature 160 K is reached by tuning 
the lattice misfit strain at fixed doping [s^. 

The calculated band structure for an ordered H 3 S lat¬ 
tice as a function of pressure clearly shows multiple Lif- 
shitz transitions for appearing of new Fermi surface spots 
in the pressure range showing high Tc superconductivity, 
which together with quantum hydrogen zero point mo¬ 
tion puts the system beyond the Migdal approximation. 

New Fermi surface spots appear at the F point at 130 
GPa pressure where the onset of the high critical tem¬ 
perature appears. It is possible that the appearing new 
Fermi surface spots drive the negative interference effect 
in the exchange interaction between multiple gaps 
contributing to the suppression of the critical temper¬ 
ature 3 . This is supported by the isotope coefficient 
which diverges at 130 GPa reaching a value of 1.5 Q, 
see Fig. [TTl well beyond the predictions of single band 
Eliashberg theory. The divergence of the isotope coeffi¬ 
cient observed here has been already observed in cuprates 
providing a clear experimental indication for a 
unconventional superconductivity near a Lifshitz transi¬ 
tion 0. 

Increasing the pressure to 210 GPa a van Hove sin¬ 
gularity crosses the chemical potential giving a Lifshitz 
transition for opening a neck. Moreover the vhs remains 


near the chemical potential within the energy range of 
the energy cutoff for the pairing interaction over the 
full pressure range between 210 and 260 GPa. We show 
that the quantum zero point hydrogen fluctuations in a 
double well Q typical of hydrogen bond and involving 
the Tu 2 phonon stretching mode, has strong effect on the 
electronic states near the Fermi level. The quantum hy¬ 
drogen zero point motion, induces fluctuations of the 600 
meV of the energy position of the vHs. The zero point 
amplitude of the S-H stretching mode, involving the Tu 2 
phonon, stabilizes the Im3m structure in the pressure 
range 130-180 GPa and induces large fluctuations of the 
small Fermi surface pockets at F . In conclusion we have 
shown the presence of large displacement amplitudes of 
ZPM. Single phonon waves can be disturbed by lattice 
quantum zero point disorder [45l| but superconductivity 
seems to resist to perturbations from ZPM in H 3 S. On 
the other hand, it is also seen that the DOS peak at 
Ep seems to pass through Ep with the large zero point 
motion. Finally more work is needed to investigate 
the variation of the Fermi level Ep in different Fermi 
surfaces with different H isotopes which change the zero 
point motion amplitude. 
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